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Abstract 

An effective method to obtain exact analytical solutions of equations describing the co- 
herent dynamics of multilevel systems is presented. The method is based on the usage of 
orthogonal polynomials, integral transforms and their discrete analogues. All the obtained 
solutions are expressed by way of special or elementary functions. 
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^ : 1 INTRODUCTION 

, Orthogonal polynomials and functions are used as the complete bases of the solutions of differential 

and difference equations. In quantum physics they are of primary importance in the construction 
of wave functions which are orthogonal by definition. In present paper orthogonal polynomials and 
functions are used to describe excitation of quantum multilevel systems in high-power laser field. 

The Schrodinger equation for a multilevel system being exciting in a monochromatic laser field 
looks like 

- ida n (t)/dt = f n+ i exp (-ie n+1 t) a n+ i (t) + f n exp (is n t) a n _i(t), n = 0,1,2,... (1) 

where f n is a dipole moment function for the transitions n— l<-m (that is fi n —i,n — /nA*o,i)j 
e n is a frequency detuning. Probability amplitudes a n (t) allow to determine levels populations 
Pn{t) = |a„(i)| 2 , an average number of absorbed photons < n >= J2n=T n Pn(t) etc. 
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To obtain the solution of (|l|) we introduce an analytical method based on orthogonal func- 
tions p n (z). 

2 FROM ORTHOGONAL POLYNOMIALS TO QUAN- 
TUM SYSTEMS 

We seek the solution of ([!]) in form 

a n (t) = ( B a(x)^P-^l exp [it (r8(x) + s n )\ dx, (2) 

J A "0 «n 

where d n is the squared norm and s n determines the frequency detuning; s n — s n — s n —\. 

At t = the solution (||) passes into the orthogonality relation for the functions p n {x). Hence 
we obtain the initial condition for Eq. (Q): a n (t) = 5 n ,o- The functions p n (x) should satisfy the 
recurrence relation 

(r0(x) + S n ) — = ./„ — h in+i— -, ■ (3) 

Initially we select the sequence of orthogonal functions which determines the quantum system 
to be investigated. The takeoff criterion is the relevant recurrence relation which allows to find 
both the constant r and the functions /„ and s n . After that we obtain the analytical solution of 
([!]) using the integral transform like (||). 

The application of analytical method (|l[^) makes possible describing dynamics of manifold 
quantum models including: A^-level and infinite-level, equidistant and non-equidistant, excited 
resonantly and nonresonantly, closed and nonclosed, excited by monochromatic field and laser 
pulse. Gegenbauer, Jacobi, Pollaczek, Laguerre, Meixner, Krawtchouk, Hahn polynomials; Szego- 
Jacobi polynomials built using the Szego method; Christoffel-Legendre orthogonal polynomials 
built by using the Christoffcl formula; orthogonal Gauss hypergeometric functions are used by us 
to obtain the exact solutions for these systems. 

This quantum systems manifold includes the systems with various characteristics, first of all 
with qualitatively different dipole moment functions /„. In particular Jacobi, Legendre and Cheby- 
shev of the kind polynomials lead to the quantum systems with f n < 1, Jacobi, Gegenbauer and 
Chebyshev of the 2 nt ^ kind polynomials — to the systems with 1 < /„ < y/n, Pollaczek, Laguerre 
and Meixner polynomials give rise to two sets of the systems with y/n < f n < n and f n > n. It is 
known that Hermite polynomials correspond to the harmonic oscillator (/„ = y/n). 

Thus the existing diversity of orthogonal polynomials makes possible to obtain the solutions of 
([!]) describing the dynamics of various quantum systems with different properties. Each orthogonal 
polynomials sequence (OPS) can be put in unique correspondence with some quantum system. 
Below we shall discuss some examples. 
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3 JACOBI POLYNOMIALS 



Let's consider the quantum systems which have been built by means of Jacobi polynomials. The 
dipole moment function following from a recurrence relation (||) looks like 



fn 



2r 



2n + a + f3 
a + (3 + 2 



n (n + a) (n + f$) (n + a + f3) 



(2n + a + /3-l) (2n + a + (3 + 1) 



a + /3 + 3 
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a > -1, /3 > 1. 



(4) 
(5) 



Thus the dipole moments of higher transitions tend to the constant value r/2. Energy levels of the 
Jacobi systems are nonequidistant but the frequency detuning decreases with n approximately as 
n~ 3 , i.e. the highest levels are practically equidistant and they are being excited resonantly. 

Jacobi OPS contains two parameters a and (3 which determine the energy levels positions. 
Therefore we obtain two-parameter continual family of multilevel quantum systems. Each pair of 
(a, 0) corresponds to the specific Jacobi quantum system which is described by dipole moment f n 
and frequency detuning e„ functions. 

The solution of ([l]) for the Jacobi systems dynamics is expressed through the Kummer hyper- 
geometric functions: 



(a+l)„ (/3+l) r 



(a + /3 + 2) 2 (n + a + [3 + l) n n\ 



(6) 



Mr+ Sn ) iFi ( y7l + a + i-2n + a + 13 + 2; -2irt) . 



At a = (3 Jacobi polynomials P( a < a '(x) pass into Gegenbauer ones CW(x), A = a + 1/2. 
Thus Kummer function in (|^) is reduced to a cylindrical function J n+ \(rt). The case a = (3 = 1 
(Chebyshev polynomials of the 2 n< ^ kind) gives so called equal- Rabi system f n = 1, for a = (3 = 
(Chebyshev polynomials of the 1 st kind) we obtain the system with f\ = l,/„ = 1/V2, n > 2. The 
levels of both Gegenbauer and Chebyshev systems are equidistant, the excitation is resonant. 

An interesting situation arises in a case a = — (3 when 



/„ = [3 (n 2 - a 2 ) I (1 - a 2 ) (An 2 l)] 
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a [3/ (l-a 2 )] 1/2 * 



n,0- 



Here the frequency detuning s n — — So^n,ij the levels are equidistant everywhere except for the 
first transition and only this transition is nonresonant. In particular at a — —f3 — ±1/2 we obtain 
/„ = 1, e n = —S n i and the expressions for levels populations p n (t) — |a„(i)| 2 reshape an especially 
simple form p n (t) = J 2 (2t) + J 2 +1 (2t). 
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4 KRAWTCHOUK POLYNOMIALS 



Using Krawtchouk polynomials of discrete variable j2J we can describe the nonresonant excitation 
of the quantum systems with finite number of equidistant levels and dipolc moment function 



fn 



n(N-n + l) 



N 



n = Q,N, N = 1,2,. 



(7) 



The analytical solution of (Q) is expressed || by way of elementary functions: 



a n (t) 




[CO 








Pn(t) 






= 










<n> 


= Ny(t) 



)y n (t) 



{[l-(l + ^M*)] 1/a +iMt)] 1/a } W n exp\ite(n-N/2)}; (8) 



iN- 



where 



y(t) 



sin 2 {[(l + a) /N] /2 t} 



l + <7 



Ne 2 /A. 



(9) 
(10) 

(11) 



The Krawtchouk systems are suitable exactly solvable models for vibrational excitation of a 
molecule in high-power infrared laser fields. The solution describes an excitation with arbitrary 
detuning of laser frequency from the system transitions frequency. The case p - 1 j 2 corresponds to 
the resonant excitation. Krawtchouk systems family includes both two-level and three-level systems 
and harmonic oscillator as special and limit cases (N = 1, N = 2 and N — > 00 correspondingly). 

Thus the Krawtchouk polynomials give a possibility to build one-parametric family of quantum 
multilevel systems and to obtain exact analytical solutions of the equations describing the excitation 
dynamics. The Krawtchouk systems family is an extensive class of multilevel models for various 
dynamical processes in spectroscopy, nonlinear optics and the other fields. 

The method under consideration gives also the possibility of solving the equations which are 
more complicated then (|l]). Let the quantum system consist of finite number of degenerate levels. 
The equation describing the resonant excitation of such a system is 

.da n , m (t) 



dt 



= n [/„+ (t) + f n a n -i 

(*)] + 



(12) 



,m+l (*) + fn+19 m&n-\-l,in— 1 (*) + fn9 m+l&n— 1,771+1 771^71 — 1,771 — 1 

(*)] 

i / i i 

where Q and fi' are Rabi frequencies, /„ = [n(N — n + 1)/N] and g m — [m(M — m + 1)/M\ 12 

are dipole moment functions for the interlevel and intersublevel transitions correspondingly, n = 

(0) 



0, N, m = 0, M; M, N = 1, 2, . . . Here we assume a n<m 
An — ±1, Am = 0, ±1 as the selection rules. 



S n ,o5 m ,o as an initial condition and 



The solution of (12) is expressed via Krawtchouk polynomials and elementary functions: 

'AT 



2 M- 



Km/ 



V2 m 

Y,€^(j,M)r) S m«T j cos N -«T j 

j=0 ^ 3 



(13) 
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where r«- 



Selecting other orthogonal polynomials we can obtain the solutions of Eq. ( |12| ) for different 
types of /„ and g m . 



5 CHRISTOFFEL-LEGENDRE POLYNOMIALS 

Let's consider the solution of (Q) using Christoffel-Legendre orthogonal polynomials 

b 



Pn{x) 



[P n ( X ) ~ K~\b)P n+1 {x)\ , K n (b) = P n+1 (b)/P n (b) 



(14) 



built with Christoffel formula j^] using Legendre polynomials P n (x). The polynomials p n {x) are 
orthogonal on an interval [—1, 1] with respect to the weight 



<j(x) = (b — x)/b, 



b> 1 



and squared norm 



d n = { [2/ (n + 1)] b/n n {b)}' 2 
They satisfy the recurrence relation ([5j) (9(x) = x) where 

fn = r(n/(2n + l))d n -i/d n , 
r = 3dx/d , 

n+1 (h\ 71 + 2 fh\ 
2n + 1 2n + 3 



(15) 



(16) 



(17) 



The application of Christoffel-Legendre polynomials in (||) gives Q the analytical solution 

'■K\h d Q exp(is n t) f . _ x -, 

[4+0,5 (^J-JK,, (b) Jn+1,5 {r t)\ 



a n (t,b)=t n [-\ ^ 



(18) 



.27 <f„ (rt)V 

of (0). This solution describes the excitation of one-parameter family of multilevel Christoffel- 
Legendre systems with nonequidistant energy levels, i.e. e n ^ const. 



6 LEGENDRE FUNCTIONS OF THE FIRST KIND 

Another quantum system can be built by using the functions 

p x n (x) = (-iJTOi + lJP^-"^) , 

/a— n 

being a special case of Legendre function of the 1 st kind Pjf (x) . The functions p^ (x) satisfy the 
recurrence relation (|J) where 

V2 



(19) 



fn 



n(2A + n-l) (A + 1) 



2(A + n—l)(X + n) 



0, 



(20) 



r = iy/2\\ + l\, 6(x) 
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and 



d n = |r(l - 2A - n)T(n + 1) /(A + n) 



A < -1, < 7i < 



(21) 



For Pn(x) is correct 



,. n r(-2A-7i + l)n! 



r (-A + i) 



A < — 1, < ra < |A| (22) 



where are Gegenbauer functions B of a complex argument z 



C X a {z) = 



T{a + 2X) 



2 F t 



T(a + l)r(2A) 
Functions are orthogonal on the discrete interval: 

N—l 



—a, a + 2A 
A + 1/2 



1-2 



J2<r iN) (xi)pi(xi)Pn(xi)=Srn,nd 2 n , JV > 1 1 — A| , A < -1 
1=0 

with respect to the weight function 
JCn 



(23) 



(24) 



(N) ( \ ( 2 A r(l - X)d N ^i \ , 2 a I A , . d A , . 



xf-1 



(25) 



where are the roots of a polynomial of degree iV, C^(zi) — 0,xi = z\ / y/ zf — 1 and /Cjv is the 
higher coefficient of this polynomial. Thus the solution of (0) for the system of N equidistant levels 
in a resonant field (at Jn = and e„ = 0) looks like 

N ~ 1 '.a^n\* x, 



(*) = E <^ W (^) f Pm(li) ) 



exp 



itr xi 



(26) 



7 FROM QUANTUM SYSTEMS TO ORTHOGONAL POLY- 
NOMIALS 

Previously we initially selected orthogonal functions to obtain the solution of the equation (|l|) for 
qualitatively and quantitatively different quantum systems. Having done so we get the compact 
solutions which are expressed via known special or even elementary functions. When OPS contains 
the parameters the solution describes the dynamics of continual or discrete family of quantum 
systems. From the other side the imperfection of this method is that it is not always possible to 
obtain exact analytical solutions for a target system with the given dipole moments function /„. 

Let's consider the quantum system with an arbitrary dipole moments function being restricted 
with the case of resonant excitation and a finite number of equidistant levels. Now (0) reduces to 
the system of N + 1 uniform linear differential equations 



ida (t)/dt = /iai(i); 



-ida„(t)/dt = f n+1 a n+1 (t) + /„a n -i(*); n = l,N-l 
-i da N (t)/dt = f N a,N-i(t) 



(27) 



(i 



(some special cases of f n have been studied in ) ■ 
Let us seek the solution of (127D in the form 



N 



l n(t) = ^ akVn ( Afe ) eXP ( iA kt) (28) 



fc=0 



where p n , A k and a k are the functions to be determined. Substituting fl28|) into (27) we come to 
a three-term recurrence relation for p n : 



fn+lPn+l (Afe) + fnPn-l (Afe) = A k p n (Afe) , 



P-i = 0, po = const (29) 



Therefore p„ are the polynomials orthogonal with respect to variable Afe (generally on a nonuniform 
grid). 

This truncated orthogonal polynomials p n (Afe) differ from the common polynomials listed above 
which satisfy the same dipole moment function f n with an interval of n change. Let's discover the 
basic characteristics of these polynomials. 

It is known that the recurrence relation like (||) determines an orthogonal polynomial up to an 
arbitrary constant factor. The comparison of (|J) and ( p9| ) shows that the polynomials p n (Afe) are 
orthonormalized. Therefore it is natural to choose the normalization po (Afe) = 1. 

Suppose that initially all the particles are concentrated on the lower level, in other words 
a n (0) = S n> o. Then it follows from ( p8| ) that 

N 

^ VkPn (Afe) p Q (Afe) = 5 n; Q. (30) 

fc=0 

Thus the integration constants a k are the weight coefficients of polynomials p n with respect to the 
variable Afe. 



To determine Afe we shall use the last equation of (27); it gives 



JnVn-i (Afe) = K kPN (Afe) . (31) 



Comparing ( |3l| ) and (29) it is possible to conclude that Afe are the roots of the polynomial 
p N+1 (Afe). 

Calculation of the polynomials p n (Afe) by the recurrence relation ( p9| ) and determination of 
weight coefficients a k from the system ([jo]) are possible for any values of N. At the same time 
( |30| ) is reduced eventually to the solving in radicals of the polynomial equations of the degree 
N+l. Generally it is possible only for A^+l < 4. However in the case of a resonant excitation the 
polynomials p n contain only even or odd degrees of an argument (depending on n parity). Thus 
the roots of the polynomials up to the ninth degree can be always found analytically. However 
using some computer algebra system (e.g. Mathematica) it is possible to obtain the roots for the 
polynomials of higher degrees for some simple /„. In any case it is possible to find these roots 
numerically with an arbitrary precision. 



7 



The method allows to obtain the solutions describing the dynamics of the quantum systems 
with an arbitrary dipole moment dependence on level number. When the solution of (Q) with the 
same /„ but for infinite number of levels can be received under the formula (|2|) by selecting the 
relevant orthogonal polynomials it is interesting to compare the truncated polynomials with the 
common ones. From the comparison of (§) and Eq. (||) it follows that 

p n (A k )=pl c HA k /r) (32) 

where pffl are common polynomials, r = 60^0/^1 f^]j d n is a squared norm, bo is defined from 
the recurrence relations for unnormalized polynomials: p n +i(x) + b n p n _i(x) = (rx + s n )p n (x). In 
particular if Rj^ +1 ' is the kfi 1 root of p^ +1 then A/. = rR^ +1 \ 

Methods ([!]-§) and @ can be generalized to the system excitation in an arbitrary laser field, 
to the account of non-adjacent dipole transitions and for the systems of degenerate levels. 
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